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Abstract. The outer layers of a neutron star are supposed to be formed of a 
solid Coulomb lattice of neutron rich nuclei. At densities above neutron drip density 
(about one thousandth of nuclear saturation density), this lattice is immersed in a 
neutron fluid. Bragg scattering of those dripped neutrons by the nuclei which has been 
usually neglected is investigated, within a simple mean field model with Bloch type 
boundary conditions. The main purpose of this work is to provide some estimates 
for the entrainment coefficients, as required for hydrodynamical two fluid simulations 
of neutron star crust which relate the momentum of one fluid to the particle 

currents of the other two fluids jS]. The implications for the equilibrium neutron star 
crust structure are also briefly discussed. 



1 Introduction 

In a typical neutron star (with mass about one and a half that of the sun), a solid phase 
of neutron rich nuclei, hence forming the crust, is expected to occur from the surface 
down to about ~ 1 km in depth where the pressure and density are so high that the 
nuclei may be strongly deformed and adopt "exotic" shapes Inside the crust, above 
the neutron drip density Pdrip — 10^^ g.cm"'^ , free neutrons are found to coexist 
with the nuclei. Below the crust, the nuclei merge into a uniform mixture of nucleons 
and perhaps of some other mesons. 

Bragg scattering of the dripped neutrons by the nuclei, which is analogous to the 
scattering of valence electrons by ions in ordinary solids, has been usually neglected 
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whereas it may be of quantitative importance for evaluating the processes involved in 
phenomena such as pulsar glitches. This issue has been pointed out by a few authors 
jS] [ZI and band effects, namely the opening of band gaps in the single particle energy 
spectrum, have been actually estimated in a simplified one dimensional case |Hj. 

In a preceeding work |2] we have introduced a macroscopic effective mass 
relevant for hydrodynamical simulations, such that the momentum of the neutron su- 
perfluid is given by the mean neutron velocity times this effective mass in the crust rest 
frame. In a (super)fluid mixture, in general the momentum of one species is a linear 
combination of the particle currents of the other components. This is known as the 
entrainment effect [H] . We have shown how to obtain these entrainment coefficients in a 
two fluid neutron star model from the effective mass 7?!^ . Our analysis is complemen- 
tary to the two fluid description of neutron star core, based on the Fermi liquid theory 
(see Borumand et al JU] and references therein), whose relativistic generalisation has 
been recently performed by Comer et al jTTj via a (7 — mean field model. 

Particle transport in the inner crust is induced by displacements of the whole super- 
fluid and it can thus be seen [H] that it is therefore not very sensitive to the superfluid 
energy gap unlike thermal properties such as the specific heat. The main effect of the 
pairing interaction is actually to lead to a smearing of the sharp neutron Fermi surface. 
Since at the densities of interest, the pairing gap is typically much smaller than the 
Fermi energy, corrections of the effective mass due to pairing correlations are expected 
to be small. We shall therefore neglect as a first approximation pairing interaction (as 
in the works of Borumand et al ^0] and Comer et al ^T] ) . 

Since the pionnering Hartree-Fock calculations of Negele and Vautherin [T^, more 
or less arbitrary boundary conditions have usually been applied. In particular the so 
called 'Wigner-Seitz approximation", whereby nuclei are treated as if isolated, has 
been widely used. Whereas this approximation would be appropriate for describing 
nucleons that are clustered around lattice nodes, it is more questionable for the dripped 
neutrons which are delocalized over the crust. This means that one has necessarily 
to consider the whole medium and apply global boundary conditions. In the case 
of a periodic lattice of nuclei, as a result of the well known Floquet-Bloch theorem 
[TB] . the problem can still be reduced to the so called Wigner-Seitz cell (not to be 
confused with the W-S sphere) namely some "elementary" polyhedron surrounding a 
nuclei whose shape is completely determined by the lattice structure, supplemented with 
Bloch type boundary conditions. The single particle energy spectrum is composed of a 
set of analytic functions £'Q,{k} (in the following we shall use brackets for functional 
dependence in order to avoid confusions with mere multiplication) of the momentum k 
(related to the translational symmetry), each sheet in momentum space being specified 
by a band index a (associated with rotational symmetry). 

Recently Magierski et al ^ have suggested that shell effects arising from "unbound" 
neutrons may be important for neutron star crust structure and have carried out Skyrme 
Hartree-Fock calculations in a cubic box with periodic boundary conditions. However 
such a choice of boundary conditions does not properly account for Bragg scattering 
since the ensuing single particle states are only those associated with vanishing Bloch 
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momentum. Besides the Wigner-Seitz cell, which is taken as a cube there, is cubic only 
in a simple cubic lattice. 

In order to prepare the way to more realistic calculations, we shall briefly review 
band theory and some key solid state physics technics which may not be familiar to 
the reader. We shall then discuss some results for dripped neutrons in the inner crust, 
described within a mean field model. 

2 Inner layer of the crust 
2.1 Mean field model 

In the absence of any previous calculation of neutron band structure except for a simple 
ID model jH], even at the simplest level of approximation, we shall adopt the single 
particle model of Oyamatsu and Yamada jS], to estimate the effective mass of the 
dripped neutrons in the bottom layers of the inner crust. Previously, we only considered 
slab shaped and rod like nuclei for simplicity The effective mass was found to be 
very close to the bare mass. The reason was that the region we were focusing on was 
nearly homogeneous. In the present work, we shall extend the same analysis for lower 
densities where the crust is assumed to be composed of spherical nuclei arranged in a 
body centered cubic lattice. Simple cubic and face centered cubic structures will also 
be considered for comparison. We shall also report results for the spherical and cylinder 
hole phases near the crust-core transition layer. 

In the model suggested by Oyamatsu et al 0, both bound and dripped neutrons 

are assumed to be described as independent particles moving in a "background" phe- 

nomenological mean field. The neutron single particle states thus obey the following 

Schrodinger equation 

-^Aif + Vif = Eif. (1) 

We further assume that nuclei are located at the nodes of an infinite crystal lattice. 
The single particle periodic potential V is thus constructed from the potential U of 
Oyamatsu and Yamada given in the Wigner-Seitz sphere by the following procedure: 

V{r} = j:U{r-T}, (2) 

T 

where the sum goes over all lattice sites, i.e. T = ^"'^a ! where are integers 

and e^j basis lattice vectors. The potential U was defined by Oyamatsu and Yamada 
as (equation (2.5) of reference [H]) 

U{r} = (^-i- j ' / dVt/o{r'}e-l-^'l^/'^^ , (3) 
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where V is & potential energy density functional of homogeneous nuclear matter. The 
local neutron and proton densities were obtained from a zero temperature Thomas- 
Fermi calculation in the W-S approximation. The energy density functional was parametrised 
such as to reproduce the properties of terrestrial nuclei on the (3 stability line and the 
equation of state of Friedman and Pandharipande for both pure neutron matter and 
symmetric nuclear matter. The parameter K (as well as the spin-orbit coupling param- 
eters, see section m below) was adjusted so as to give reasonable values for the single 
particle energies of ^^^Pb (for further details, see [SHE])- We have used the parameter 
sets of model I from reference jH]. We have neglected the spin-orbit coupling introduced 
by Oyamatsu and Yamada since it is smaller by an order of magnitude compared to the 
central potential as shown on figure ^ In order to take into account Bragg scattering, 
single particle states have to satisfy the Floquet-Bloch theorem 



As a consequence, the Schrodinger equation (P) can be solved within any primitive cell 
[T^ , among which the Wigner-Seitz one possesses the full symmetry of the lattice, with 
boundary conditions given by (0) • The Wigner-Seitz (W-S) cell is defined as the set of 
points that are closer to a given lattice node than to any other. Exemples are shown for 
cubic crystals on figured Except for simple cubic lattice, the W-S cell is a complicated 
polyhedron. In the work of Oyamatsu and Yamada, the lattice spacing was defined as 
Q'Oy — "^cif 5 where V^^jj is the volume of the W-S sphere. In the present work, we 
have defined the lattice spacing a for cubic lattices (cube length) such that the volume 
of the W-S polyhedron is equal to that of the sphere. 

2.2 Single particle energy spectrum 

For each given k , equation with boundary conditions (0) will admit a discrete set of 
eigenvalues £^q, (indexed by a Greek letter) referred as bands. Single particle energies 
will be ordered such that Si < S2 < ••• which guarantees that for a given band, 
say a, the energy ^Q,{k} will be continuous (besides it is analytic except for high 
symmetry points). Unlike the W-S approximation, the single particle energy spectrum 
in the Bloch case has a much more detailed structure. 

It can be shown that the single particle energy is periodic in a "reciprocal" lattice 

[T!^ whose nodes are given by K = eLi ^al" ' where ia are integers and the dual 

basis is defined by the dot products 



(/.k{r + T} = ei^V{r}. 



(5) 



r ■eb = 27T5%. 
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Figure 1: Single particle central potential (left) and spin-orbit coupling potential (right) 
in the W-S sphere for the lowest baryon density considered = 0.03 fm from 
model I of reference jH]. 




Figure 2: Wigner-Seitz cell of cubic lattices (from left to right) simple cubic, body 
centered cubic, face centered cubic with conventional labelling of high symmetry points 
and lines [22] • 



Figure 3: Nuclear configurations taken from reference jH^: (a) spherical nuclei in a 
body centered cubic lattice, (b) "spaghettis" in 2D hexagonal lattice, (c) "lasagna" , (d) 
"antispaghetti" in a 2D hexagonal lattice and (e) "bubbles" in a body centered cubic 
lattice. 

This means in particular that equation (^J) has only to be solved for each momentum 
k , within a domain known as the first Brillouin zone (BZ), i.e. the W-S cell of 
the reciprocal lattice. By considering rotational symmetry, it can be shown that all the 
relevant information about the spectrum is contained within some "irreducible" domain 
of the BZ . The solutions of Magierski et al [U] are only those associated with the 
center of the first Brillouin zone, namely k = 0, provided the lattice is simple cubic 
(for body centered or face centered cubic structures, their solutions are not of the Bloch 
type). 

We have solved the Schrodinger equation with Bloch type boundary conditions 
dni) variationally by a finite plane wave expansion of the single particle wavef unctions, 
namely 



i:^k{K}e'*+'^>-- 



(7) 



with an energy cutoff £^cutoff such that 
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(8) 



The nuclear structures that we have considered are illustrated on figure El 
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2.3 Effective mass 



The zero temperature many body ground state is merely obtained by filling all the 
single particle states whose energy £^Q,{k} is lower than the Fermi energy /i . For 
each total neutron density TT-n, the Fermi energy /i is thus determined via the integral 
over the first Brillouin zone (BZ) 
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where we have explicitely introduced the Fermi distribution by the Heaviside unit step 
distribution = 1 if X > 1 and zero otherwise, allowing each momentum state 

to be g -fold degenerate (here g = 2 to account for the spin degeneracy). In particular, 
the ground state will be characterized by the Fermi surface bounding the occupied part 
of momentum space. In general the Fermi surface will consist of disconnected pieces, 
each piece being given by ^Q,{k} = /i. 

Introducing the dripped or free neutron density n , the effective mass (from which 

the entrainment coefficients can be obtained is defined P by the simple formula 

TL 

= ^, (10) 
where the mobility scalar JC is given by (summation over energy bands is implicit) 

and the group velocity is defined by the usual expression 

V = iVk^. (12) 
h 

This mobility scalar is related to the electric conductivity (7 by the formula (7 = e^T JC 
in which e and T are the electric charge per particle and some characteristic relaxation 

time respectively. In terms of the effective mass, this formula reduces to the well known 

expression [13] (T = ne^r/m^. 

In the following sections, we have defined the dripped neutron states to be those 

states whose energy is positive £^Q,{k} > (the energy origin being taken as the 

largest possible value of the potential). 



2.4 Wigner-Seitz approximation 

The W-S cell of body or face centered (unlike simple) cubic lattice is a nearly spherical 
polyhedron, this is the origin of the so called W-S approximation according to 
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which the W-S cell is taken as a sphere and more or less arbitrary boundary conditions 
are applied on the single particle wave functions. For instance, the prescription given 
by Negele and Vautherin ^2] requires that on the sphere radius, radial wave functions 
with even orbital angular momentum quantum number / vanish and derivatives of the 
radial wave functions with odd / vanish. Besides the density needs to be averaged in 
the vicinity of the cell edge in order to prevent the occurrence of unphysical density 
fiuctutations. Within the current framework, the Floquet-Bloch theorem (0) ensures 
that the local total neutron density defined by 

^n{r} = -^E/ d3fc|^.k{r}|M/i-^a{k}} (13) 
possesses the lattice symmetry, the single particle wavefunctions being normalized as 

^/ dhMv}\' = l. (14) 

y „ "'cell 

coll 

This means in particular that the density gradient VT7n{T} will vanish whenever 
the crystal is left invariant under space inversion, which means that the density will 
be nearly constant (to first order) inside nuclei. In cases for which r belongs to 
a symmetry plane, the density gradient will have no components orthogonal to this 
plane. This means that the density profile will be essentially fiat (to first order) around 
the cell boundary, along directions perpendicular to the W-S faces whenever the face 
is parallele to any symmetry plane. This is true for the simple cubic and face centered 
cubic W-S cells. This is also true for the body centered cubic W-S cell excluding the 
hexagonal faces. 



3 Numerical methods 
3.1 Brillouin zone integrations 

Brillouin zone integrations are involved in many solid state physics calculations, as for 
instance in the particle density (fT3|) . Analytic expressions for the integrand are usually 
unknown (except for a few academic cases), therefore one has to rely on numerical 
schemes based on discrete summations. Since the computation of energy bands is 
usually the most time consuming part of solid state calculations, several techniques 
have been developed in order to compute these integrals with the smallest number of 
terms. One of the most successful methods, which was pioneered by Baldereshi and 
subsequently improved by Chadi-Cohen jT^j, is based on a weighted sum over a set of 
special symmetry points. The main idea relies on the fact that a function /{k} which 
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has the full symmetry of the reciprocal lattice, can be expanded over symmetrized plane 
waves (in the following we shall consider only cubic structures, for which there is no 
glide plane or screw symmetry axis [32]) 

/{k} = E/.p^Ee^''^^•^ (15) 

s \i \ V 

where "P is any rotation of the lattice and \P\ their number (lattice vectors are 
ordered as [T^J > l^sal whenever S\ > S2 )• The integral of this function over 
the first Brillouin zone is rewritten as a sum over a set of unequivalent points kj with 
weight iy{kj} : 

1 r ^ 

fo = ^J /{k}d'^ = E^{kj}/{kj} + i?„ (16) 
E^{kj} = l (17) 

where the residual is given by 

~ N 1 

i?. = - E /. E ^{kj}p^ E e^'''^^-'^ • (18) 

By choosing the points kj such that all plane waves are exactely zero for all S 
up to some largest Smax ? the residual will contain only terms with lattice vectors 
|T| > IT^^^^I for which the contribution will be vanishingly small provided the 
function /{k} is sufficiently smooth. Eventually the integral is approximated by 

1 . N 

-— / /{k}d3A::^E^{kj}/{kj}. (19) 

l^BZ j = l 

The determination of these special points kj and their associated weights w;{kj} was 
rather complicated in the first scheme involving a recursion process. Several authors 
have provided formulae (see ^Sj and references therein). 

The convergence of the discrete sum with respect to the number of special points 

is exponential for smooth integrands and an error of less than one per cent can be 

attained within a few points only. Unfortunately the presence of a sharp Fermi surface 

(especially in any metallic solid) introduces discontinuities in the integrands, such as in 

equation (|T3|) , which spoil the convergence. This method is thereby most satisfactory 
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Number of special points 

Figure 4: Relative error in the computation of integrals (j2T| with respect to the number 
of points 



for insulating materials and semiconductors. Nevertheless, accurate results can still be 
obtained with a small number of special points by "smearing" the Fermi surface P^ . 
Since the volume of the first Brillouin zone and of the Wigner-Seitz cell are related by 

Ve. ^ ^. (20) 

the number of special points for a given precision is smaller for larger cells. We have 
tested the special point method by calculating the integral 

k'd'k, (21) 

which is equal to 1/4, 3/8, 19/32 in units of (27r/a)2 for simple cubic, body cen- 
tered cubic and face centered cubic lattice respectively. The convergence is illustrated 
on figure |3] We have also shown on figure |S1 how the smearing procedure of Methfessel 
et al [T^ could improve the convergence in the extreme case of one single special point, 
by computing the total neutron density as a function of the Fermi energy n defined 
in equation (jH)) in the Shockley "empty" lattice test [201 which we have 



10 




Figure 5: Density as a function of the Fermi energy for the "empty" lattice model 
(body centered cubic lattice, lattice parameter=30 fm), with one single special point 
and with different smearing parameters {N is the number of approximants and W is 
the smearing width, see reference ^Hl for the details) 



3.2 Fermi energy 

The determination of the Fermi energy is numerically expensive as it requires the cal- 
culation of all occupied single particle states (unlike the effective mass which depends 
only on the highest occupied states on the Fermi surface). The computional cost is 
significant since the number of levels is of the order of several hundreds up to nearly 
one thousand. Nonetheless a useful estimate of the Fermi energy with an accuracy of 
a few percents can still be easily obtained as follows. First, among all occupied states, 
one should distinguish between core states whose wavefunction is localised in the neigh- 
borhood of nuclei and valence or conduction states whose wavefunction extends over all 
space. The core states are easily obtained by finding the bound states of one isolated 
nucleus. Having found the number of such states, let's say A^core ; equation (jH)) can be 
written as 



Tin = 7/73 VezA^corc + 7/73 E / ^'k ^fl - £:«{k}}l?{i'«{k} - ^-eore}, 

(23) 

where Score is the energy of the highest core state. The determination of the Fermi 
energy is essentially a state counting procedure which becomes less and less sensitive to 
the detailed structure of the energy spectrum as the number of particles is increased. 
Consequently assuming that the valence states can be roughly described as non inter- 
acting particles, the previous equation thus leads to the following estimate: 
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where from (pUj) we have TT-core = 5'^core/"^cell- 

For instance, for the baryon density = 0.03 fm , we found for the core states: 
three s states, two p states, two d states, one f state and one g state. The total number 
of core states is thus equal to A^core — 35 and the Fermi energy is approximately 
given by /i ~ 18.02 MeV, which is about 5% larger than the value /i = 17.16 
MeV, obtained by direct integrations of via special points. The value of the Fermi 
energy was also found to be completely unsensitive to the lattice structure, which is 
due to the fact that equation (jH)) does not explicitely depend on the precise values of 
energy bands but only on their number in a given energy range. 



3.3 Fermi surface integrations 

Calculations of the effective mass (fTUI) require a fine mesh sampling of the Fermi surface. 
One of the most accurate methods for performing Fermi surface integrations is based 
on the Gilat-Raubenheimer (GR) scheme [25 . The integration over the Fermi surface 
is first reduced by symmetry to the irreducible wedge of the first Brillouin zone which 
is partitionned into microcells. The integral is thereby discretized as a sum over the 
microcells with suitable weight factors to account for the fact that microcells may overfill 
each wedge. The Fermi surface is then approximated by a plane inside the microcells 
which it intersects, so that the integration within each cell can be calculated analytically 

We have implemented the original linear extrapolation method with cubic cells j^H] ■ 
The energy is linearly extrapolated from the values of the energy and velocity at the 
center of the cell. This scheme therefore requires the energy bands as well as their gradi- 
ents, whose evaluation with plane waves is straightforward via the Hellmann-Feynman 
theorem (2^, namely 

v = — + E^^k K 25 
m \^ TTi 

with the normalisation of the wavefunction given by 

El^k(K)P = l. (26) 

K 

The GR method seems to be not so popular as the tetrahedron method in solid state 
physics, in which the irreducible domain is divided into tetrahedra and the energy is 
linearly interpolated from its values at the four corners of each tetrahedron |25] |26j . 
One of the main reasons is that the calculation of the velocity may be computationally 
expensive with more elaborate basis functions than plane waves. Nevertheless the 
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linear extrapolation of the energy within the GR method, which is a first order Taylor 
expansion around the cell center, is a better approximation than an interpolation. Since 
in the present case the integrand in equation flll|) involves the energy gradient, the GR 
method seems much more appropriate. Besides interpolation faces the band crossing 
problem: the ordering of the energy bands with increasing energy ensures that the 
energy is continuous for a given band but does not prevent discontinuities in the energy 
gradient from occuring. Unlike extrapolation, interpolation will yield an incorrectly 
small gradient resulting in a loss of convergence (see [2H] [21] for a discussion about 
systematic errors). 

In order to test this scheme, we have employed the "cubium" toy model, i.e. a 
simple cubic crystal with one single s band whose momentum dispersion in cartesian 
coordinates is given by ^{k} = — (cos + COS ky + COS A;^). Analytic expressions 
for the density of single particle states at energy S , defined by 



are available in terms of the elliptic integral of the first kind ^U]- As can be seen 
on figure ^{g = 2 ) the GR method yields results in very good agreement with the 
analytic expressions, in particular in the vicinity of the critical points associated with 
V = Q (producing the "kinks" in the density of states). 

4 Results 

We have considered the three types of cubic lattices, namely simple cubic, face centered 
cubic and body centered cubic for spherical nuclei and bubbles {i.e. periodic lattice 
of "holes" containing dripped neutrons in uniform nuclear matter). The Fermi surface 
integrations have been carried out with of the order of ~ 10*^ cubic cells, for which 
the relative deviation was a few percents as shown on figures [7| and |H1 

The effective mass shown on figurelTUlhas been found to be greatly enhanced due 
to Bragg scattering, reaching ~ 15 772 at a total baryon density TT-b — 0.03 fm 
Deeper inside the crust, the renormalization is smaller, tending to ~ m at 
very high density as the single particle potential tends to a constant. We have also 
shown on figure ^2 the transverse effective neutron mass in the "exotic" phases PP 
including the "antispaghetti" layer, where non spherical nuclei occur (only the velocity 
components respectively perpendicular to the slabs or to the rods, for the "lasagna" 
or "(anti) spaghetti" phases are included in equation (|TT|l for the mobility scalar ; the 
other velocity components have the same expressions as those of a noninteracting gas 
and therefore the corresponding effective mass coincide with the ordinary mass). The 
dependence on the lattice structure is rather weak. Numerical results are summarized 
in appendix ^ We have also carried out the calculations with the estimate of the 
Fermi energy suggested in section 13.21 While this is usually a good approximation for 




(27) 
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Figure 6: Density of states of the cubium model computed with the GR method and 
compared with the analytic expressions given by Morita et al [201 
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bulk quantities such as the conduction neutron density, the errors are much larger for 
quantities depending only on states at the Fermi level. In particular, the resulting 
effective mass is about 20% lower at baryon density = 0.03 fm , having value 
?7i^/ TTl — 14 , where the relative error in the conduction neutron density is less than 
on percent. 

In order to understand the origin of this effective mass enhancement, we have com- 
pared the Fermi surface area Sp with its expression obtained from the assumption 
that the dripped neutrons form an ideal Fermi gas of density n , which is given by 

5^, = 47r^(67rV^)2/3. (28) 

We have found that the Fermi surface area is strongly reduced in the outer layers 
compared to that of the sphere, whereas it tends to its uniform expression ()28|1 near 
the crust-core interface as illustrated on figure El It is first to be remarked that the 
"enclosed" Fermi volume depends only on the density, as is easily seen on equation 
© and therefore does not depend on the particular shape of the Fermi surface. It is 
well known that among all the closed surfaces of a given volume, the sphere has the 
minimum surface area. In the present case however, the opening of numerous band gaps 
in the single particle energy spectrum leads to a Fermi surface composed of disjoint open 
pieces. By band gaps, we mean avoided crossings of enery band sheets in k space. 
This is the reason why the Fermi surface area is found to be smaller than that of the 
sphere in the low density layers of the crust where the gaps take their largest values. 

It is quite instructive to evaluate the mean Fermi velocity Vp defined by 

Vp = ^fvdSp. (29) 

Op F 

It is readily seen that it is expressible, in terms of the Fermi surface area S and 

^ ' gas 

Fermi velocity v of the non interacting neutron 



S m 

F gas 



(30) 



gas 



For instance, at the lowest density we considered, namely TT.^ = 0.03 fm , the mean 
Fermi velocity is about one third the Fermi velocity of the gas. The ratio increases to one 
half at TT-b = 0.055 fm and tends to one as the density approaches the crust-core 
transition density. This is a clear and intuitive manifestation of Bragg scattering. 

The spin-orbit coupling, as we have seen, is very small compared to the central 
potential in the bottom layers we have considered (see figure Q). However such a term 
would probably not be negligible in the lower density regions of neutron star crust. The 
presence of this term in the Schrodinger equation preserves the lattice symmetry 
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which imphes that single particles wave functions can still be expressed as Bloch states. 
However the spin-orbit coupling breaks the spin symmetry and therefore raise some 
degeneracy. As a result, this will open new gaps in the energy spectrum, in the vicinity 
of which the group velocity will be vanishingly small. Consequently, it can be seen from 
(II Ij) and (fTn|) that the resulting effective masses would be expected to be increased by 
the spin-orbit coupling since the density would not be much affected. This conclusion 
is confirmed by a numerical calculation in the "spaghetti" phase. The computations 
are quite involved since the size of the secular (complex) matrices is doubled. We have 
taken into account the spin-orbit coupling given by Oyamatsu and Yamada (equation 
(2.12) of reference [S]), which is given in the W-S approximate cell (a cylinder in this 
case directed along the Z axis) by an expression of the form 




where r is the distance to the cylinder axis, the angular momentum component 
along the cylinder axis and the diagonal Pauli spin matrix. The parameters Ai and 
A2 were adjusted so as to give the correct sequence of single particles states of ^'^^Pb . 

We have used the values of model I from reference [S] namely Ai = 175.8 MeV.fm^ 

and A2 = 16.39 MeV.fm^ . We have carried out the calculation at the total baryon 
density = 0.06 fm for the hexagonal lattice. Since the densities varies smoothly, 
the spin-orbit coupling has a small effect, yet the opening of band gaps can be clearly 
observed as shown on figure E| The transverse effective mass is found to be slightly 
larger with than without spin-orbit coupling as expected, respectively T\l^ j ni = 1.37 

and mj^ /m = 1.35 while the Fermi surface area is smaller, respectively S^/ S^_^^ = 
0.88 compared to S^/ S ^ = 0.89 . 



5 Conclusion 

We have considered the band effects in neutron star crust induced by Bragg scattering of 
dripped neutrons upon the solid lattice of nuclei. This is the analog in a nuclear context, 
of a well known phenomenon in ordinary solid state physics whereby valence electrons 
are being scattered by ions leading notably to an understanding of semiconductor band 
gaps. Based on a mean field model with Bloch type boundary conditions, we have 
computed the effective dripped neutron macroscopic mass necessary for evaluating 
the entrainment coefficients in hydrodynamical descriptions of neutron stars ^1 12] . We 
found that this effective mass is increased by several hundred percents compared to the 
ordinary neutron mass in the middle layers of the inner crust but becomes negligible 
near the crust-core boundary where nuclei merge into a uniform mixture of neutrons. 
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Figure 7: Convergence of the effective neutron mass ?7l^/ m with respect to the number 

of GR cells, for baryon density ?2b ~ 0.055 fm with nuclei in a body centered cubic 
lattice 
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Figure 8: Convergence of the Fermi surface area S / S with respect to the number of 

GR cells, for baryon density TT-b = 0.055 fm with nuclei in a body centered cubic 
lattice 
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Figure 9: Neutron Fermi surface area S^/S for configurations with spherical nuclei 
(filled symbols) and bubbles (unfilled symbols) in cubic lattices 
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Figure 11: Transverse effective neutron mass m^/m for "exotic" configurations 




Figure 12: First Brillouin zone and irreducible domain of an hexagonal lattice with the 
conventional labeling of high symmetry points and lines. 
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M r KM 

Figure 13: Neutron band structure, with (thin hues) and without (thick hues) along 
high symmetry directions shown on figure [T^ for an hexagonal lattice of spaghetti (rod) 
like nuclei at baryon density Tli, = 0.055 fm 




Figure 14: Zoom of the region in the box of figure E| Band splittings induced by 
spin-orbit coupling can be clearly seen. 
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electrons and protons (and perhaps of some other mesons). We have also found that 
the lattice structure plays a minor role. However the differences are significant between 
3D crystals and the ID or 2D liquid crystal configurations that were previously studied 
within the same model (respectively slab shaped and rod like nuclei). 

The present results also suggest that the widely used Wigner-Seitz approximation, 
whereby Bragg scattering is neglected, may be questionable for the study of neutron 
star crust. In particular, it is to be remarked that the very large values of the effective 
mass that we found in the middle layers tend to indicate that dripped neutron 
shell (actually band) effects may have significant consequences for the neutron star crust 
equilibrium structure compared to that obtained in the W-S approximation. It has been 
recently shown that the composition of the inner crust of neutron star is quite sensitive 
to the proton shell effects [22] • The unbound neutron shell effects have been investigated 
by Magierski and Heenen [Bj who suggested that such quantum effects play an important 
role in the determination of the neutron star crust equilibrium structure since the 
energy difference between different nuclear shapes, lattice structures and compositions 
is typically small. Their analysis was based on a semiclassical approach. They have 
recently carried out Skyrme Hartree-Fock calculations [HI] but with ordinary periodic 
boundary conditions. The implementation of Bloch type boundary conditions within 
more realistic models, especially taking into account neutron pairing as outlined recently 
0, should therefore be further investigated. 
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APPENDIX 



A Numerical results 

Effective dripped neutron mass for splierical nuclei and bubbles in a body centered 
cubic lattice: 



nb(fm ^) 


77n(fm ^) 




m^/m 


0.03 


0.029 


0.94 


15.4 


0.055 


0.053 


0.95 


5.5 


0.0855 


0.0826 


0.98 


1.07 



Transverse effective dripped neutron mass for "pasta" phases: "spaghettis" (hexagonal 
lattice), "lasagna" and cylinder holes (hexagonal lattice): 



nb(fm 3) 


^n(fm ^) 


n/n^ 


m^/m 


0.06 


0.0581 


0.9553 


1.347 


0.065 


0.063 


0.9578 


1.242 


0.07 


0.0678 


0.9606 


1.174 


0.074 


0.0716 


0.9633 


1.124 


0.076 


0.0735 


0.9634 


1.031 


0.0775 


0.0749 


0.9646 


1.029 


0.08 


0.0773 


0.9666 


1.025 


0.082 


0.0792 


0.9687 


1.020 


0.083 


0.0802 


0.9734 


1.010 


0.084 


0.0811 


0.9746 


1.003 


0.085 


0.0821 


0.9762 


1.001 
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